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Abstract 


A series of flight tests to measure the wake of a Lockheed C-I30 
airplane and the accompanying atmospheric state have been conducted. 

A specially instrumented North American Rockwell OV-IO airplane was 
used to measure the wake and atmospheric conditions. An integrated 
database has been compiled for wake characterization and validation of 
wake vortex computational models. This paper describes the wake- 
measurement flight-data reduction process. 

Introduction 

Several of today’s major airports are operating their capacity limit, leading to an increase in airport 
congestion and delays. As more and more airplanes are placed into the terminal area the probability of 
encountering wake turbulence is increased. A wake vortex upset is most hazardous for aircraft near the 
ground during landing and takeoff. The degree of upset mainly depends on the relative sizes of the vortex 
generating and vortex encountering airplanes, and the extent of wake decay. The rate of wake decay is 
highly dependent on the atmospheric state. 

Numerous laboratory experiments and analytical models have been developed to study and simulate 
wake vortex flow physics. However, there is very little full-scale data available for comparison with and 
validation of the experimental and computational results. Many of the previous wake-measurement flight 
tests have omitted the atmospheric state data, which has a direct influence on wake flow physics. The 
National Aeronautics and Space Administration (NASA) conducted a series of flight tests to develop a 
wake-measurement data set with the accompanying atmospheric state information. This data set has been 
compiled into a database that can be used by wake vortex researchers to compare with experimental and 
computational results. Details of the database are provided in reference 1. This report describes the 
wake-measurement flight-data reduction process. 

Symbols 

a 

b C130 

b a 

b v 

h APq 

^AA/p, ’ h APqp-> w * n g' boom position error pressure correction due to sideslip Dias terms, psi 
b ^ flank angle position error calibration bias, deg 

C a 5-hole probe angle of attack calibration curve coefficient 


speed of sound, fps 
wing span of C- 1 30, ft 

angle of attack position error calibration bias, deg 
sideslip position error calibration bias, deg 
position error pressure correction bias, psi 

nose-boom position error pressure correction due to sideslip bias, psi 


t * t r i f 
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5-hole probe sideslip calibration curve coefficient 


east velocity integration correction factor 

north velocity integration correction factor 

vertical velocity integration correction factor 

dynamic pressure calibration coefficient 

static pressure calibration coefficient 

5-hole probe angle of attack coefficient 

5-hole probe sideslip coefficient 

distance east of wake origin, ft 

east distance computed from GPS, ft 

east distance computed from integrated inertial velocity, ft 

wake origin east offset distance, ft 

distance north of wake origin, ft 

north distance computed from GPS, ft 

north distance computed from integrated inertial velocity, ft 

wake origin north offset distance, ft 

wake drift distance, ft 

vertical distance below wake origin, ft 

vertical distance computed from GPS, ft 

vertical distance computed from integrated inertial velocity, ft 

wake origin vertical offset distance, ft 

longitudinal wake drift distance in wake axis system, ft 

lateral wake drift distance in wake axis system, ft 

airplane inertial velocity component east, fps 

gravitational acceleration, 32.2 ft/s 2 
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n gmd 
h pass 

HI 

i 


K APq 


AP % 

k APc^ ’ k APc^2 


*1 

K 


altitude, ft 

altitude at the ground, ft 
altitude of a pass, ft 
GPS geometric altitude, meters 
array index variable 
array index variable 

linear slope of angle of attack position error calibration 
linear slope of sideslip position error calibration 
linear slope of position error pressure correction 

linear slope of nose-boom position error pressure correction due to sideslip 
linear slope of wing-boom position error pressure correction due to sideslip 
linear slope of flank angle position error calibration 
temperature recovery factor (0.995) 


Lat 

latitude, deg 

Lat a 

latitude of point a , deg 

Lat b 

latitude of point b , deg 

U ' C r,f 

latitude of C- 1 30 at t re j , deg 

Lat N 

latitude north, deg 

lMt 0 

latitude of wake origin location 

^OVref 

latitude of OV- 1 0 at i re j , deg 

LatoV starl 

latitude of OV-IO at t s(arr deg 

ny 

u y stop 

latitude of OV- 1 0 at t slop , deg 

Lng 

longitude, deg 

L »g a 

longitude of point a, deg 

Ln Zb 

longitude of point b , deg 
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u ' g c„ ; 

longitude of C- 1 30 at t re f , deg 

L "io 

longitude of wake origin location, deg 

u ' s ov re/ 

longitude of OV-10 at t re f , deg 


longitude of OV-10 at t siarr deg 

l "8 0 V , 
w stop 

longitude of OV-10 at t stop , deg 

LngW 

longitude west, deg 

M 

Mach number 

Mag var 

magnetic variation 

N 

airplane inertial velocity component north, fps 

Outfile 

Wake_Vel output file name 

Pm 

roll rate from mechanical gyro, deg/sec 

PDOP 

dilution of precision parameter 

p <> 

standard sea-level static pressure, 21 16.22 psf 

P ref 

reference pressure, 1000 millibars 

Ps 

static pressure, psi 

Ps c 

probe calibration corrected static pressure, psi 

Ps o 

static pressure corrected for zero bias, psi 

p i 

5-hole probe center port pressure, psi 

P al 

5-hole probe bottom port pressure, psi 

P a 2 

5-hole probe top port pressure, psi 

P P. 

5-hole probe right port pressure, psi 

P P2 

5-hole probe left port pressure, psi 

Q 

dynamic pressure, psi 

4c 

probe calibration corrected dynamic pressure, ] 

Qm 

pitch rate from mechanical gyro, deg/sec 
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<!» 

R 

RMS 

Rug 

R "8err 

R "Smin 

SV 


'GPS 


1 pass 


start 


1 stop 


f ref 


abs 


l be 


T 

oo 

H 


dynamic pressure corrected for zero bias, psi 
yaw rate from mechanical gyro, deg/sec 
gas constant, 1716.5 ft lb/(slug °R) 
root mean squared position error, meters 
range from one point to another, n. mi. 
range error, n. mi. 
minimum range error value, n. mi. 
number of navigation satellites received 
time, sec 

global positioning system time, sec 

time when encountered wake was generated, sec 

time of pass, sec 

event start time, sec 

event stop time, sec 

time of vortex encounter, sec 

absolute total air temperature, °R 

body axis to Earth axis transformation matrix 

Earth axis to wake axis transformation matrix 

total air temperature, °C 

standard sea-level temperature, 518.67 °R 

potential temperature, °R 

free-stream air temperature, °R 

body-axis longitudinal airspeed component, fps 

Earth-axis longitudinal airspeed component, fps 

body-axis lateral airspeed component, fps 
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Vrtx 

Vt 

V '!oc 

V wnd 

w b 

w e 

W P 


w 


N 


Wnd Dir 


Wnd, 


Wnd 


N 


Wnd 


Spd 


Wnd, 


Wild y 


w? 


x, v, 


x b'>b' z b 


* y o ' ^s) 

U W U W U W 


ov t 


stan 


-ov 4 


stop 


a 


Earth-axis lateral airspeed component, fps 
vortex identifier (L for left, R for right) 
true airspeed, fps 
local true airspeed, fps 
wind speed, knots 

body-axis vertical airspeed component, fps 
Earth-axis vertical airspeed component, fps 
Earth-axis easterly wind component, fps 
Earth-axis northerly wind component, fps 
wind direction, degrees magnetic 
easterly wind component, fps 
northerly wind component, fps 
wind speed, knots 

longitudinal wind component in the wake axis system, fps 
lateral wind component in the wake axis system, fps 
Earth-axis vertical wind component, fps 
airplane axis coordinates, ft 

airplane axis coordinates of the body-axis origin, ft 
wake axis coordinates of the body-axis origin, ft 
geometric altitude, ft 

geometric altitude of wake origin location, ft 
geometric altitude of OV-IO at f starf , ft 
geometric altitude of OV-IO at t stop , ft 
airplane inertial velocity component up, fps 
angle of attack, deg 


7 



a Vn 

P 

P/ 

P Vn 


Y 

T1 

P 

P grnd 


P mid 
Po 

e 


¥ 

'F 


w 

wnd 
A Pq 


Ap q$ 


A R 


a/2 


AFp/2 


AP 


3a/ 


^3(3/ 

At 


angle of attack corrected for airplane rotation rates, deg 
wing tip alpha vane measurement, deg 
sideslip angle, deg 

sideslip angle corrected for airplane rotation rates, deg 

wing tip beta vane measurement, deg 

specific heat ratio for air, 1 .40 

flank angle, deg 

dummy variable 

air density, slug/fr 

air density at the ground, slug/ft 1 

air density at middle altitude, slug/ft 3 

standard sea-level air density, 0.0023769 slug/ft 

pitch angle, deg 

roll angle, deg 

yaw angle, deg 

heading, deg 

heading of wake axis system, deg 

wind direction, deg 

position-error pressure correction, psi 

position-error pressure correction due to sideslip, psi 

pressure difference between 5-hole probe bottom and top ports, psi 

pressure difference between 5-hole probe right and left ports, psi 

pressure difference between 5-hole probe center and bottom ports, psi 

pressure difference between 5-hole probe center and right ports, psi 

wake age, sec 
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Subscripts 


abs 

absolute value 

avg 

average value 

Bias 

zero bias correction 

Cl 30 

C-130 value 

e 

Earth axis system 

Fit 

from least-squares linear fit 

grnd 

value at the ground 

l 

left boom 

mag 

degrees magnetic 

n 

nose boom 

pass 

value during a pass 

probe 

probe measured value 

r 

right boom 

ref 

reference value 

tip 

C-130 wing tip value 

w 

wake axis system 

wk 

wake value 

wnd 

wind 

5 

5-hole probe 

93 

Litton LN-93 INU value 

Superscripts 


raw, uncorrected measurement 

Abbreviations 


a/c 


aircraft 



CPT 

control position transducer 

DAS 

experimental data acquisition system 

INU 

inertial navigation unit 

GPS 

global positioning system 

LAFB 

Langley Airforce Base 

NACA 

National Advisory Committee for Aeronautics 

NASA 

National Aeronautics and Space Administration 

PCM 

pulse code modulation 

OVDRA 

OV-10 data reduction and analysis program 

SCADC 

standard central air-data computer 

UTC 

Coordinated Universal Time 

WakeJVel 

wake velocity program 

Flight Test Overview 


The objective of the wake-measurement flight tests was to develop a full-scale data set with the 
accompanying atmospheric state information. The tests used two NASA airplanes as illustrated in 
figure I. 

The NASA Wallops Flight Facility’s Lockheed C-130, shown in figure 2, was the wake generator. It 
was outfitted with wing tip smokers to mark the wake. It weighed between 105,000 and 95,000 pounds 
during the test and has a wingspan of 132 feet 7 inches. The tests were all flown with the C-130 in the 
clean (flaps up, gear up) configuration. 

The NASA Langley’s North American Rockwell OV-10A measured the wake and atmospheric 
conditions. The OV-10 was equipped with a three-boom, flow-sensor arrangement to measure the flow 
field characteristics of the wake. The booms were located on each wing tip and the right side of the nose 
as shown in figure 3. The OV-10 instrumentation is described in more detail in the next section. The 
OV-10 flew through the wake at various downstream distances measuring the wake velocities and 
position. A picture of a wake taken from the video camera mounted on the left tail is shown in figure 4. 
The smoke clearly denotes the two vortices. 

The general flight test procedure was as follows. The OV-10 would fly to the flight test area and 
begin a series of “weather” runs. These would begin with one or two “turbulence” runs, which were 2- 
minute, level, constant heading flight segments 1000 feet below the test altitude. The test altitude was 
generally 4000 or 5000 feet. This was followed by a “weather profile” run, which was a constant speed, 
constant heading climb from 1000 feet below to 500 feet above the test altitude. Another series 
turbulence runs would then follow 500 feet above the test altitude. After completing the weather runs the 
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OV-IO would rendezvous with the C-130 and begin a series of “wake measurement” runs. The C-130 
would fly a constant speed and heading at the test altitude. The OV-IO measured the wake by flying 
slower than the C-130 and making a series of wake penetrations at increasing ranges behind the C-130. 
At the conclusion of a series of wake measurement runs an additional series of weather runs were then 
conducted. 

OV-IO Instrumentation 

Significant modifications to the OV-IO were required to convert it from a stock aircraft into a research 
platform. An overview of the development of the research systems and measurement capabilities of the 
OV-IO is presented in reference 2. The following section provides a summary of the research 
instrumentation and experimental data system. 

General System 

A data acquisition system (DAS) was developed to collect and convert the analog signals from the 
experimental sensors and convert them into digital form for display, transmission, and recording. The 
system could record approximately two hours of data from 150 parameters at data rates up to 128 Hz. 
The major experimental instrumentation components are labeled in figure 3. The sensor outputs were 
signal conditioned, multiplexed and digitized by a commercial 12-bit pulse code modulation (PCM) 
subsystem that features functional flexibility, programmability, and the ability to multiplex analog and 
digital signals into a serial format. Major components of the sensor suite include: 

• time-code generator to provide an accurate time base for all flight measurements, 

• ring-laser-gyro inertial navigation unit (INU), 

• global positioning system (GPS), 

• standard central air-data computer (SCADC), 

• fast-response rate and attitude gyros, 

• two pitot-static and flow-angle-measurement systems, 

• 5-hole airspeed and flow-angle probe, 

• temperature probes, 

• dew-point sensor, 

• video/audio cameras/recorders, 

• control position transducers (CPTs). 

Except for sensors, displays, and control panels, most of the equipment for the DAS system is installed 
in the OV-IO cargo bay. A removable main instrumentation pallet, designed to fit in the aft portion of the 
cargo bay and roll out onto a cart for system maintenance, was fabricated and populated with equipment 
(Figures 5a and b). The forward part of the cargo bay housed the majority of the navigation, avionics and 
display equipment (Figure 6). 

Navigation Systems 

Three satellite based navigation systems have been added to the OV-IO. The first is a multi- 
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component system that supports the research measurements and includes a Litton LN-93 ring-laser gyro 
INU integrated with a Honeywell 3A GPS receiver and a CPU- 140/A SCADC. All the components 
communicate with the PCM system via a MIL-standard 1553 data bus through an 80486-processor/33- 
MHz personal-computer bus controller. The INU updates the GPS in the event of satellite signal loss, but 
the GPS does not update the INU position estimate. The INU also is capable of providing aircraft 
accelerations, velocities, and attitudes. Supplemental measurements of the linear accelerations and 
angular rates are provided by a separate package of fast-response rate gyros and accelerometers. The 
components of the research navigation system reside in the forward cargo bay on fixed plates and are not 
part of the removable instrument pallet. The second GPS system is an Ashtech Z12 system capable of 
storing 4 hours of position and velocity data at 1 Hz data rate. The third navigation system entails a 
Garmin GPS 100 mounted in the forward cockpit, and enables the pilot to precisely navigate and/or set up 
research maneuvers. 

Air-Data Systems 

There are two air-data systems on the OV-10. One system is the SCADC previously described, which 
is coupled to the INU/GPS system. The second system is the three-component pitot, static, and flow- 
direction measurement system described here. 

Sensors mounted on the tips of graphite-epoxy research booms facilitate air data measurements at both 
wingtips and the nose. As is customary, the booms were designed to place sensors as far ahead of the 
predicted aircraft upwash influence as possible, yet have enough lightness and stiffness (high natural 
frequency) to prevent unwanted vibration influence on the sensor measurements. Thermally-controlled 
Setra pressure transducers, located approximately at the base of each respective research boom, provide 
for enhanced accuracy in pressure measurements. Accelerometers mounted near the tips of these booms 
provide information on boom vibration. 

A standard National Advisory Committee for Aeronautics (NACA) pitot-static probe with balsa angle 
of attack and sideslip vanes, shown in figure 7, was mounted at the end of each wingtip boom 
approximately one chord length in front of the wing leading edge. The static and total pressure ports were 
connected to the pressure transducers by 13 feet of 1/4-inch inside-diameter flexible tubing. The balsa 
vanes, which have a high natural frequency of about 0.18 Hz per knot of indicated airspeed in 
incompressible flow (reference 3), have a natural frequency of approximately 22 Hz at the nominal 
OV-10 airspeed of 120 knots. 

A Rosemount Model 858 AJ 5-hole probe, shown in figure 8, was mounted at the end of the nose 
boom. The Rosemount probe provides the required pressure ports for determining flow direction, 
airspeed, and altimetry measurements at the nose. The five pressure ports at the tip of the probe and the 
ring of static pressure ports about 3.5 inches aft are connected to a manifold such that the following 
pressure differences between the ports shown in figure 8 are measured. 


AP(xI2 — Pa\ Pa2 

(psi) 


(1) 

<N 

1 

£ 

II 

£ 

(psi) 

X 

(2) 

^3a 1 ~ p 3~ p a\ 

(psi) 


(3) 
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AP 3PI = P 3- p P! 


(psi) 


(4) 


q n =P 3 -Ps n (psi) (5) 

The tip of the mounted probe is approximately 6 feet in front of the aircraft nose and approximately 2 
feet to the right of the longitudinal centerline. The probe was connected to the pressure transducers by 
approximately 17 feet of 3/1 6-inch inside-diameter flexible tubing. A small 0.067-inch inside diameter 
restrictor was inserted in each 3/16 tubing near the 5-hole probe. The restrictor size was experimentally 
selected per reference 4 to provide maximum damping to the pressure oscillations measured by the 
transducers when a 20 Hz sinusoidal pressure was applied to the probe. 

A pair of Rosemount Model 102 non-deiced total temperature probes (one for the SC ADC and one for 
the research system) and a General Eastern 1011 Aircraft Dew-Point Sensor are also integrated into the 
air-data system. The temperature probes are located on the underside of the wing approximately midway 
between each engine and wingtip. The dew-point sensor is mounted on the left side of the forward 
fuselage. 

Airborne Video/Audio Systems 

The research instrumentation includes a video/audio system consisting of three miniature lipstick-size 
Elmo video cameras and three recorders. One camera is located at the tip of the left vertical tail, pointing 
in the general direction of flight, and includes the top of the wing and fuselage in its field of view, as 
shown in figure 4. This camera provides a good, qualitative means of reviewing the highlights of a flight 
in preparation for data reduction and analysis. The other two cameras are mounted near the tip of each 
wing and are pointed vertically downward. The wing tip cameras provide a stereo recording capability 
that is post-processed to compute the spacing between the vortex smoke trails and distance from the 
OV-10. Audio channels on the video/audio system enable the pilot and flight test engineer to record 
verbal data during maneuvers. 

Flight Data Reduction Process 

The basic data flow of the flight data reduction process is shown in figure 9. From each flight four 
data products are generated directly from the OV-10. 

• Time-tagged, video recordings from the tail camera and both wing-tip cameras. 

• Flight notes from the test pilot and flight engineers of event times and conditions 

• Time-tagged data tape from the experimental data system 

• Ashtech differential GPS data 

The flight notes and video data are reviewed and used to establish a precise log of event times and test 
conditions called the “flight event” file. This information is then used to extract from the experimental 
data tape the pre and post-flight instrument calibration data and data from special calibration flight 
maneuvers. This calibration data is used to establish the instrumentation zero bias values, which are 
saved as a “flight constants” file for each flight. The event and constants files and the data from the DAS 
are input for the OV-10 Data Reduction and Analysis Program (OVDRA). This program applies the 
instrument calibration and bias corrections to the DAS data and computes the inertial referenced wind 
components. Details of the OVDRA algorithms are presented in the Calibration and Bias Corrections 
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section of this report. The output of OVDRA is a calibrated data file for each run of the test flight. These 
run data files are used to characterize the atmosphere through wind and temperature profiles and 
turbulence calculations. 

The wake velocity measurements are computed from the OVDRA output, the flight event file and the 
GPS position data from the OV-IO and the C-130. This information is input to the Wake Velocity 
Program (Wake_Vel) which computes the origin and age of the wake and translates the wind 
measurements to the wake axis system. Details of the Wake_Vel algorithms and the wake axis system 
are provided in the Wake Velocity Calculation section. The output of the Wake_Vel program is a data file 
for each wake measurement event of each run of each flight. 

The pilot notes from the C-130 are used to estimate the initial circulation strength of the wake, which 
is required for subsequent wake decay analysis. 

Calibration and Bias Corrections 


The instrument calibrations and bias corrections are applied to the DAS data through the OVDRA 
program. All the recorded DAS data parameters are listed in table 1 along with their symbol and database 
name. All the data are interpolated to the highest sampling rate of 128 Hz. The sampling rates of all 
parameters are listed in the table. Those data parameters that do not have a symbol representation in 
table 1 are not used in the data reduction algorithms. 

One of the primal^ functions of the OVDRA program is to translate the OV-IO atmospheric 
measurements from the airplane or body axis reference system, shown in figure 10 to an inertial reference 
system. This transformation is made through the yaw (if/), pitch (9), and roll (<j>) Euler angles shown in 
figure 1 1 . 

The OVDRA program reads the DAS data at each time step and applies the following calculations. 
Zero bias corrections 


First the zero bias corrections are applied to the following parameters. The bias values were 
determined from pre-flight measurements with the airplane stationary in the hanger. These corrections 
compensate for shifts produced by the DAS signal conditioning circuitry. 


• The static and dynamic pressure measurements from the left, nose and right booms, respectively. 


r \, """ 

j-l,n,r 


( 6 ) 


Qoj = Vj' -Vjtiias ( P Si) 
j—l,n,r 


(7) 


• The pressure difference measurements across the nose boom 5-hole probe. 

r 

Mai 2 = APcf/2 _A/ a/2 Bias (ps>) 


( 8 ) 
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Wpi 2 =&Ppi 2 - AP j3/2 Bias (psi) 


(9) 


AP 3a ,=APj a , -A/> ia/Bjas (psi) (10) 

AP 3 pi =APjp, -AP 3 p, ms (psi) (11) 

• The mechanical gyro measurements of roll, pitch and yaw rates, respectively. 

/ 

Pm~Pm ~ Pm Bj as (deg/ sec) (12) 

/ 

Qm-Qm ~~ ^WRias (deg/sec) (13) 

r m = r m “ 'm Bias (deg/ sec) ( 1 4) 

• The total temperature measurement. 

Tm = T m jas ( ^ (15) 

Rosemont 5-Hole Probe Calibration Correction 


The first step in the 5-hole probe calibration correction is to compute the angle-of- attack and sideslip 
coefficients 


, ^ AP o(\2 

2 + )+ 4 { AP P\ 2 + A/ al 2 ) 


(16) 


C fl = 


AF/ 


^12 


9 (aF 3 « i + AP^I )+ 1 (aP« , 2 + AP a | 2 ) 


(17) 


A fifth-order two-dimensional polynomial calibration curve that is a function of C a and C p is then 
used to compute the local angle of attack and sideslip. 


4 4 


«5 = XIC'C' C (deg) 

(=o )=o 


(18) 


4 4 


05 = X X Ca C’n -Ck. (deg) 

1=07=0 


(19) 


The local dynamic and static pressure are computed in a similar manner. 
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(psi) 


(20) 


45 


q _% 

4 4 . . 
J-XlQrC' 
i=07=0 


C 95/+7+! 


^5 = Ps 0 +<7 5 X I C& CL ■ C„ (psi) (2 1 ) 

/=0 7=0 

The values for the polynomial coefficients C fl , C/ 7 , and C s were derived from proprietary 
manufacturer’s data and independent wind tunnel tests. 

Position Error Corrections to Local Flow Angle Measurements 

The position error correction terms were derived from the airspeed calibration tests that are discussed 
in appendix A. The corrections to the local angle of attack measurements are: 


a r ~ ka r Vn f . + b a r 

(deg) 

(22) 

a / = a Vni 

(deg) 

(23) 

a n a 5 + 

(deg) 

(24) 


The vertical balsa vanes on the wing tip booms that are referred to as the sideslip or P vanes actually 
measure the flank angle rather than the true sideslip. The position corrections to the flank angle 
measurements are: 


r ” *r| r -Pv, V + b x\ r 

(deg) 

(25) 

*1/ “*rj/ ’Pvtoy + b X] t 

(deg) 

(26) 


The sideslip angles are computed from the flank angle measurements through the trigonometric 
relationship. 


P j = tan 1 (tan r? j cos aj ) (deg) 
j = r,l 

The position correction for the nose boom sideslip measurement is: 

P« =*P„ P.s + &p n ( de 8) 


(27) 


(28) 


The baseline values for the position error correction terms (i.e. the k and b terms) are provided in table 2, 
which lists all the calibration constants. The b terms are modified for each flight based on data from pre 
and post-flight calibrations and in-flight calibration maneuvers. 
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NACA Pitot Probe Calibration Corrections 


The NACA pitot probe measurements from each wing tip boom are corrected for local flow angle 
effects based on wind tunnel calibration data. The static and dynamic pressure correction coefficients, 
C s and C q respectively, are obtained from a 2-D table lookup as a function of the local angle of attack 

and sideslip. Tables 3 and 4 list the C s and C cj coefficient values, respectively, as a function of the local 

angle of attack and sideslip. The dynamic and static pressure corrections are: 


qc j 1 -c 


q j 


j = r,l 


Ps c j = Ps Oj +C Sj q cj 


(29) 


(30) 


j-r.l 

Position Error Corrections to Static and Dynamic Pressure Measurements 

The position-error pressure correction terms were determined from the airspeed calibration tests and 
are developed in Appendix A. The correction terms were found to be a function of the local dynamic 
pressure and sideslip. The pressure correction terms are computed as: 


AP( 1 j = k A P qj ■ Q Cj + h AP<ij + APc i pj (P si ) 

j = r,l 

A p cj„ = k A p (/n ■ q 5 + b APqn + APq K (psi) 
where the sideslip correction is: 


(31) 

(32) 


= Min^O, A 

AP( f$ n *|P«l + ^AP(/p w 

)■ 

l^APqp i ' 

Pr + ^>APqp\ 

> Pr 

<0 

\ k APqp2 ' 

Pr +b APqp 2 

. Pr 

>0 

\~kAPqp 2 

Pi + b APqp 2 

.Pi 

<0 

| ~^APqp t 

Pi + b APqp\ 

. Pi 

>0 


(psi) 

(psi) 

The position error pressure correction is applied to the static and dynamic pressure through: 

Ps j = Ps c j - APq j (psi) 

j = r,l 

Ps n = Ps 5 -APq n (psi) 


(33) 


(34) 


(35) 


(36) 

(37) 
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q j =qc j +AFq j (psi) 

j = r,l 

(38) 

< 1 ,, =45 +A ^4n (psi) 

(39) 


Local Mach Number Calculation 

After all the corrections are applied to the static and dynamic pressure measurements the local Mach 
number at each boom is computed from: 


y-1 


Mj = 


y-1 


±L_ 

Psi 


+ 1 


-1 


(40) 


j = n,r,l 


where the ratio of specific heats (y) for air is 1.40. The average local Mach number is then computed as: 

Mm-# =ji{M n + M r + M i) 

Free-Air Temperature and Speed of Sound Calculation 

The measured temperature is converted from degrees Celsius to absolute temperature in degrees 
Rankine. 


T ah s= §T W +32+459.67 (°R) 


(42) 


The free-stream air temperature is then computed from the local absolute temperature as: 

Tabs 


T =- 

1 oo 


1 + 


Y ' \K M 2 aYg 


^ ( R) 


(43) 


where the temperature recovery factor (K) is 0.995. The potential temperature is computed as: 


T 0 = T a 


' p ref $ 
Ps n 68.9475 


(°R) 


(44) 


where the reference pressure ( P re f ) is 1000 millibars. 


The free-stream speed of sound is computed from: 


a = 


P 0 T„ 

V Po Tq 


(fps) 


(45) 
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where: 


P 0 = 2116.22 (Ib/ft 2 ) 
p 0 =0.0023769 (slug/ft 3 ) 
r,, =518.67 (°R) 


The local true airspeed is computed from the local Mach number and the free-stream speed of sound. 


Vt locj = Af j' a (fps) 
j = n,r,l 


( 46 ) 


Airplane Angular Motion Correction to Body Axis Airspeed Components 

The body-axis airspeed components are computed and corrected for the airplane angular motions as: 

u bj = Vtlocj cos ® j cos P j G j ~~ yb ) — Qm (*- j ~ z b )] (fp s ) 

j =n,r,l 


V. bj = Vtlocj sin Pj+Jw t j ~ Zb r,n ( X J ~ Xb )1 (fps) 


(47) 


(48) 


j = n f r t l 

w bj ~ Vt i oc j sin a j cos /Jy + — ^ \q m (a y — x y> )~ p m Gy ~ y b )] (fp s ) 
j - n, r, l 


(49) 


The airplane axis x, y and z coordinates of the boom probes and body axis origin ( , Z( ? ) are listed 
in Table 5. The angular rates (Pm*<im' r m) are from the mechanical gyros and are sampled at 32 times 

per second while the velocity parameters are sampled at 128 times per second. The angular rates are 
interpolated in the missing time slots (3 points between each sample). The angular rates from the 
mechanical gyros were selected over the LN-93 INU rates due to a dithering of the INU rate signals. 

True Airspeedy Angle Of Attack and Sideslip Calculation 


The true airspeed at each boom is computed as: 

Vt j=-ffij +V lj +W bj ( f P s ) 
j = n,r,I 


(50) 


The angle of attack at each boom is computed as: 


a, . = — tan 1 
' J n 

j = n,r,l 


( u . \ 
% 

V U »J J 


(deg) 


( 51 ) 
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Calculation of Wind Components in Earth Axis System 

The airspeed components are transformed from the body axis system to the Earth axis system through 
transformation equation: 

u bj 

v h 
% 

j = n,r,l,avg 

where the transformation matrix ( T/ Je ) is defined as: 

„ coswuisineq^sintbui cosy 93 sin 0 93 cos<J> 93 

cosy 93 cos0 93 _ sin ; 3 93COS 9 4 3 +s i„ V93 sin$ 93 

. llf siny 93 sin0 93 sin0 93 siny 93 sin0 93 cos<)> 93 

V93 ■ 93 +c0Sl|/ 93 C0S<j> 93 — cosy 93 sin <|> 93 

-sin0 93 cos0 93 sin<|> 93 cos0 93 cos<|) 93 

The inertial wind components are computed as the difference between the airplane inertial velocity 
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and the measured airspeed components. 


W Nj =N-u e . (fps) 
j = n,r,l,avg 


W Ej =E-v ej (fps) 
j = n,r,l,avg 

W z = -Z - w (fps) 

j e J 

j = n,r, l,avg 

A verage Wind Speed and Direction Calculation 

The average horizontal wind components are computed as: 


w N avi; +W N/i +W Nr ) (fps) 

W E(n . g +W En + W Ef ) (fps) 


The wind speed and direction are then computed from: 


UW W = 05925 


M +W l ( knots ) 

y avt* w/vj? 


Wnd Di r = 1 80 + tan' 


f Wf A 


W N 


avfr 


(deg true) 


(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 


At the conclusion of the above calculations the DAS input data and selected computed parameters are 
recorded. Table 6 lists the recorded parameters and cross-references the equation where they are 
computed. The data for the next time step is read and the calculation process begins again at equation 6. 
This process loop is repeated until the end of the DAS file is reached. 

Wake Velocity Calculation 

The following process was used to derive the wake velocity measurements from the OVDRA output, 
the flight event information, and the GPS position data from the OV-10 and the C-130. The objective of 
the wake measurement flight tests was to obtain a “real world” data set for use in wake modeling analysis. 
The ideal data set for such analysis would be a series of instantaneous, 3-dimensional wake measurements 
within a fixed airmass, from initial wake rollup until final decay. However, currently there is no way to 
make such measurements. The OV-IO could only measure the wake velocity at 3-points in space at any 
instant in time. A wake measurement pass consists of 3 streams of wake velocity data collected at 128 
Hz, over a 5 to 10 second period. Successive wake measurements were made in different airmass regions 
rather than an ideal fixed airmass. The OV-10 and C-130 inertial positions are known from differential 
GPS measurements collected at 1 Hz. The following process was used to translate the measured wake 
velocity data to an inertial-referenced frozen axis system. 
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Wake Coordinate System 

The wake vortex axis system (X W ,Y W ), shown in figure 12, is aligned with the bearing between the 
OV-IO and C-130 at the time of the wake vortex encounter {t ref ). The origin of the coordinate system is 
the C-130 position when it generated the wake being measured (C-130 position at t a ). 

Wake Velocity Algorithm 

The wake velocity algorithm WakeJVel must translate the OVDRA wind measurements to the wake 
coordinate system. Each wake measurement event has it own wake coordinate system. As noted in the 
previous section, the origin of the coordinate system is the C-130 position (latitude, longitude, and 
altitude) when it generated the measured wake event. Since the wake translates with the airmass this 
origin location must be iteratively determined based on the known C-130 track and OV-IO measurement 
location, and the average wind component determined from measurements outside the wake. The wake 
origin location also determines the wake age. 

Filtering GPS position data. The wake velocity algorithm begins by reading and filtering the OV-IO 
and C-130 Ashtech differential GPS data files for a given flight. The read and filter operations are the 
same for both the OV-IO and C- 1 30 data. The following data parameters are read at each time step: 

tops GPS time, seconds 

SV Number of satellites received 

PDOP Dilution of Precision Parameter 

Lufyv Latitude north, degrees 

Lngw Longitude west, degrees 

HI Geometric altitude, meters 

RMS Root mean squared position error, meters 

The following criteria are used to filter poor position data from the GPS files. The data that do not meet 
the criteria are omitted. 


RMS < 1 .0 

PDOP < 4.0 (68) 

SV>4 


For the remaining data, the GPS time is converted to Coordinated Universal Time (UTC) through: 

r GPS _ * 6.0 for 1 995 flights 
r C p 5 -11.0 for 1996 flights 

1 = ' t Gp$ - 1 2.0 for 1 997 nights 1 1 ' 

t C ps - 26.7 for flight 705 (DAS time not reset) 

The longitude west is converted to longitude east, and the geometric altitude is converted to feet through: 

Z = -3.281*77/ (ft) (70) 

Begin wake velocity computation loop. The wake velocity computation loop begins with reading the 
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event data. The loop is repeated until the end of the event data is reached. The event information consists 
of: 

t start Event start time, UTC seconds 

t st0 p Event stop time, UTC seconds 

t re j Time of vortex encounter, UTC seconds 

Vrtx Vortex identifier (L for left, R for Right) corresponding to t re j 

V wru j Wind speed, knots 

^wtul Wind direction, degrees true 

Outfile WakeJVel output file name 

Find OV-IO locations corresponding to event start, stop and reference times. The first step in the 
computation loop is to interpolate the 1 Hz OV-IO GPS data to find the location corresponding to t start » 
t stG p , and t re f . The algorithm will not interpolate between points with data dropouts. If the time interval 

between a pair of GPS data points is greater than 1.2 seconds, then the data between those points did not 
meet the dropout check criteria. An error message stating that the GPS data is invalid over the 
interpolation region is written and the event output file is closed. The algorithm returns to the beginning 
of the loop to read the next event data. 

Determine when C-130 generated the wake encountered. This is the beginning of a search-loop to 
determine when the C-130 generated the wake encountered. The search begins at the C-130 location 
corresponding to t re j and searches backward to find the time (t a ) when the wake encounter position 

(corrected for the wind drift) is closest to the C-130 track. A flowchart of the search algorithm is shown 
in figure 13. 

t o = ref ’ Lat OV re j ’ Ul SOV ref ’ Vrtx ' wnd > V wnd ’ C " 1 30 GPS data ( 7 1 ) 

The algorithm, as depicted, finds the nearest C-130 data point corresponding to the wake generation 
point. The algorithm has a second iteration loop (not depicted for clarity) which interpolates between 
GPS data points to find the wake origin time to the nearest tenth of a second and corresponding location 
(, Lat o ,Lng o ,Z 0 ). The search algorithm also uses several functions relating range and bearing between 

two latitude and longitude points which are not explicitly expressed in figure 13. These functions are 
given below. 

The range from (Lat (n Lng a ) to {Lat ^ y Lng ij) is computed as: 

where ^ = sin Lat a sin Lat ^ +cos Lat a cos Lat cos {Utg^ -Lng a ) 

The bearing from {Lat a , Lng a ) to {Lat [ 7 , Lng ) is computed as: 
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iy/ _ 1 80 
JC 


y + tan 1 




v ’ / 


(deg) 


where = 


sinLafi -sin La/ cos 
^ « ou 

t Rng 

cos Lat, t sin -~- 
u 6U 


(73) 


The latitude and longitude of a point at a given range ( Rng in nautical miles) and bearing (T) from an 
initial location (Lat a , Lng a ) are computed from: 


Lot = Lat a + cos T (deg) 


L/ig - Lng a +- 


/?/ig sin T 


/ . ^ 
Lar + Lat n 


60 cos 


[l-l-W^Rng j 


(deg) 


(74) 

(75) 


Compute wake age. At the conclusion of the search for t a the wake age is computed from: 


^ ? ref ^ q ( 76 ) 

Find C-130 location corresponding to event start and stop times. The C-130 GPS data is interpolated 
to find its location at t start - At and t stop - At . As with the OV-IO GPS interpolation, the interpolation 

interval is checked for validity. If the time interval between interpolation data points is greater than 1.2 
seconds, then the GPS data did not meet the dropout filter criteria over that region. An error message 
stating that the GPS data is invalid is written and the event output file is closed. The algorithm returns to 
the beginning of the loop to read the next event data. 

Find C-130 location at the reference time. A linear extraploation of the C-130 ground track from 
r start “ Af to t sta p -At is used to compute its location at t re f . The linear projection approach is used 

rather than interpolating the GPS data because at the end of the runs the C-130 often turns or enters a 
circular pattern until the OV-IO completes its wake measurements. Interpolating the GPS data for the 
C-130 position in such a situation will result in the wake axis system being misaligned. 

Least squares linear fits of the C-130 GPS latitude and longitude as functions of time are computed 
from t start - At to t stop - At . The projected latitude (Latc ref ) and longitude (Luge ) corresponding to 

t re f are then computed from the linear fit. 

Compute vortex initial location at the reference time. The vortex initial location at t re f is assumed to 
be the left or right wing tip location of the C- 1 30. The wing tip location at t re j is computed as: 

Lot Up . Lng Up = fi La, C llrf , Ln SC„. eJ - 'Vfip >\ h CI30 ) ( 77 ) 


using equations 74 and 75, respectively; where 
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( 78 ) 


\^CJ30 +90 , Vrtx = R (right vortex) 

{i P j 4^7 jo -90 , Vrtx - L (left vortex) 

and 

^CUO = fV^’ Cl 30 {t start ~ Al l L "S C 7 3o{ ( start _A? ) Lal C U<)k stop “H Ln gcuokstop “^)J C 79 ) 

from equation 73. 

Compute heading of wake axis system. The bearing from the OV-10 location at t re f and the C-130 tip 
location at t re f is the heading of the wake axis system. The wake axis heading is obtained as: 

4\ v = / (L at tip , Dig tip , Lat ov ref > btg OV ref J ( 80 ) 

from equation 73. 

Compute Earth to Wake axis transformation matrix. The Earth to wake axis system transformation 
matrix elements are computed as: 

r cosT w sin X V W 0^ 

T w = -sin 4^, cos x ¥ w 0 (81) 

0 0 1 

Compute OV-10 INS velocity integration correction factor The OV-IO position at each 128 Hz data 
record is determined by integrating the LN-93 INS velocity components over the duration of the event. 
An integration correction factor is then calculated that forces the position calculations to match the 1 Hz 
Ashtech GPS position data. 

The north, east and vertical distance traveled during the event are first computed by integrating the 
LN-93 velocities from t start to t stop . 


l stop 

d N lNS = X i 


-*i - 1 

) (ft) 

(82) 

1 start 





{ stop 

d Z,NS ~ ^ "2^‘ 

-7-lx £ ; 


) (ft) 

(83) 

1 start 





Hoi 

Jh 

ii 

N 

^3 

7 - 7_i \z 

i-Zi - 1 

) (ft) 

(84) 


start 


The range and bearing from the event starting location to the stop location are then computed from 
equations 72 and 73, respectively. 
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Rng, V - f\L*'OV slar , . ^80V mrl - ^'OV 5Mp . ^OV w „„ ) 


(85) 


The north, east and vertical distances based on the OV-IO Ashtech GPS data are computed as: 


dN GPS = 6076. \Rng- cos 4* 

(ft) 

(86) 

dE GPS = 6076. 1 ■ Rng ■sin'V 

(ft) 

(87) 

dz c,PS = z ovw, ^ ov stop 

(ft) 

(88) 


The north, east and vertical integration correction factors are then computed. 


Cor : 


d Jars d Jins 

t stop start 


j = N,E,Z 


(89) 


Compute the event origin offset distances. The origin of the wake axis system is the location of the 
C-130 at tg (equation 71). The range and bearing from the event starting location to the C-130 initial 

location are computed from equations 72 and 73, respectively. 

Rng ,4 ' = / 1 Latov sl(irl • L "8 OV starl ,Lat a , Lng 0 j (90) 


The origin-offset distances are then computed as: 

d N = 6076. 1 ■ Rng cos T (ft) 

1 O 


(91) 


d F = 6076. 1 Rng sin T (ft) 


(92) 


= Z 


OV start 


(ft) 


(93) 


Compute OV-IO position and wind velocities in the wake axis system. The OV-IO position in the 
wake axis system is computed at each data record by integrating the INS velocities with the integration 
correction factors and origin offset. The measured wind velocities are also transformed to the wake axis 
system. 

The INS velocity components are integrated with the correction factors and origin offsets for 

t start — 0 — 1 stop '■ 


d Nj ~ d N 0 + Eor N {>i f start ) + X \ (fj 1 j-\\^ j ^ j- 1) ^ 


(94) 


J 1 start 


d Ej = d E 0 +Cor E ('i - 1 start )+ X j H ^ j ~ E j-\ ) ( ft > 


(95) 


J-' start 
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( 96 ) 


d z i ~ d Z Q +Cor Z i'i start ) + X _ 4 ( r j ~ l j-\ K Z j ~ Z j-\ ) ( ft) 

j~ l start 

where i and i start are the data array indices corresponding to and t star(9 respectively. 

The distance components are then transformed to the wake axis system yielding the coordinates of the 
OV-IO body axis origin in the wake axis system. 


f \ 

-v 


M 

y °w 

~ T 

l ew 

dE 

\ U ) 

i 

K d Z l 


(ft) 


The boom coordinates in the Earth axis system are computed. 


) 

, f \ 

Xj e 


X j x h 

y Je 

II 

K 

yj -yt 

i J 

t 

\ J / 


(ft) 

[ z J~ Zh J 

j = r, /, n 

where the components of the matrix are obtained from the OVDRA file at each time step. 
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The boom coordinates and wind measurements are then transformed to the wake axis system. 


X 


( \ 


r \ 

X Jw 


x o w 


X ie 

y Jw 

— 

y o w 

+ T w 

y Je 

i 7 h- j 

i 

k J 

i 

yj? J 


j = r, /, n 


(ft) 


X 

“Jw 


( \ 
Wnj 

V J w 

- T 

l ew 

W E j 

J 

i 

W Z , 

\ J Ji 


j = r, l, n 


(fps) 


( 99 ) 


(100) 


Wake position correction for wind drift. This correction is made in an effort to obtain an instantaneous 
sample of the wake velocity flow-field. Since a wake measurement pass requires a number of seconds to 
complete, the measurement positions must be corrected to account for the drift of the wake over the 
measurement period. 

First, the average wind, which is one of the inputs to the Wake_Vel algorithm, is translated to the 
wake axis system. 
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where 



( 104 ) 


Write wake velocities and positions to output file. The wake velocities and positions at each boom are 
written to the output file Outfile. 

Closing Remarks 

This paper has described and documented the flight-data reduction process for a series of flight tests 
that were conducted to measure airplane wake and atmospheric state information. A specially 
instrumented North American Rockwell OV-IO airplane was used to measure the wake and atmospheric 
conditions. A Lockheed C-130 airplane was used to generate the wake. An integrated database has been 
compiled for wake characterization and validation of wake vortex computational models. 
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Table 1 . Continued. 
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Table 2. Calibration coefficient values. 


Symbol 

Default 

Value 

Units 

Description 

ka r 

0.8223 

- 

Linear slope of right vane angle of attack calibration 

V 

1 .0073 

- 

Linear slope of right flank angle vane calibration 

{ 

0.8363 

' 

Linear slope of left vane angie of attack calibration 


0.9999 

- 

Linear slope of left flank angle vane calibration 

h a r 

-1.7568 

deg 

Calibration bias for right angle of attack vane 

V 

1.4417 

deg 

Calibration bias for right flank angle vane measurement 

b «, 

-0.8308 

deg 

Calibration bias for left angle of attack vane 

% 

-0.1695 

deg 

Calibration bias for left flank angle vane measurement 


0.8002 

- 

Linear slope of the nose boom angle of attack calibration 


0.9183 

- 

Linear slope of the nose boom side-slip calibration 

b «n 

0.4420 

deg 

Bias term of the nose boom angle of attack calibration 

% 

0.5686 

deg 

Bias term of the nose boom side-slip calibration 

^APqn 



0.08849 j 

| 

Slope of nose boom position error correction curve for zero side-slip 

^APqn \ 

-0.00752 i 

J j 

psi 

Y-intercept of nose boom position error correction curve for zero side-slip 

^A Pqr ] 

j 0.014903 1 

! | 

| - ; 

Slope of right boom position error correction curve for zero side-slip 

APqr 

-0.00236 ] 

psi 

Y-intercept of right boom position error correction curve for zero side-slip 

k-APql j 

0.018634 

| - j 

Slope of left boom position error correction curve for zero side-slip 

^APql | 

-6.00362 

psi 

Y-intercept of left boom position error correction curve for zero side-slip 

,, | 

-0.00294 

j 

] 

psi/deg 

Slope of nose boom position error correction due to side-slip 

b 

W Pi j 

0.000661 ; 

psi/deg 

Slope of wing boom position error correction due to side-slip 

*AA/p 2 j 

0.000670 1 

i i 

psi/deg 

Slope of wing boom position error correction due to side-slip 


i 0.010261 | 

i J 

psi 

Y-intercept of nose boom position error correction curve due to side-slip 

&AP *pi j 

! o 

s -i 

psi 

Y-intercept of wing boom position error correction curve due to side-slip 

b &P % 2 \ 

0 

psi 

Y-intercept of wing boom position error correction curve due to side-slip 
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Table 3. NACA pitot probe static pressure calibration coefficient ( C s ) values. 



101, deg 

a, deg 
-40 

0 

5 

10 

15 

20 

25 

30 

35 

40 

0.2976 

0.3071 

0.3245 

0.3403 

0.4263 

0.5335 

0.6680 

0.6153 

0.5247 

-38 

0.2673 

0.2761 

0.2922 

0.3088 

0.3771 

0.4839 

0.6534 

0.6481 

0.5440 

-36 

0.2482 

0.2562 

0.2709 

0.2880 

0.3407 

0.4473 

0.6194 

0.6651 

0.5635 

0.5833 

-34 

0.2325 

0.2405 

0.2554 

0.2729 

0.3147 

0.4121 

0.5817 

0.6648 

-32 

0.2158 

0.2248 

0.2411 

0.2594 

0.2963 

0.3717 

0.5456 

0.6387 

0.6027 

-30 

0.1960 

0.2067 

0.2253 

0.2453 

0.2815 

0.3360 

0.481 1 

0.6077 

0.6181 

-28 

0.1763 

0.1885 

0.2127 

0.2348 

0.2723 

0.3203 

0.4270 

0.5716 

0.6262 

-26 

0.1518 

0.1648 

0.1981 

0.2221 

0.2617 

0.3104 

0.3892 

0.5455 

0.6291 

-24 

0.1308 

0.1439 

0.1785 

0.2040 

0.2458 

0.2976 

0.3600 

0.5149 

0.6252 

-22 

0.1144 

0.1272 

0.1564 

0.1843 

0.2277 

0.2827 

0.3366 

0.4814 

0.6151 

-20 

0.0942 

0.1065 

0.1372 

0.1697 

0.2151 

0 2729 

0.3235 

0.4020 

0.6006 

-18 

0.0553 

0.0673 

0.1 190 

0.1599 

0.2084 

0.2685 

0.3182 

0.3646 

0.5839 

-16 

0.0366 

0.0482 

0.0975 

0.1499 

0.2014 

0.2634 

0.3127 

0.3588 

0.5676 

-14 

0.0264 

0.0.378 

0.0746 

0.1348 

0.1878 

0.2513 

0.3025 

0.3493 

0.5538 

-12 

0.0195 

0.0307 

0.0619 

6.1258 

0.1791 

0.2437 

0.3005 

0.3465 

0.5424 

-10 

0.0143 

0.0252 

0.0539 

0.1 190 

0.1725 

0.2375 

0.3008 

0.3479 

0.5332 

-8 

0.0097 

0.0204 

0.0495 

0.1150 

0.1688 

0.2337 

0.3012 

0.351 1 

0.5258 

-6 

0.0062 

0.0167 

0.0477 

0.1130 

0.1675 

0.2318 



0.3010 

0.3544 

0.5197 

-4 

0.0033 

0.0138 

0.0467 

0.1115 

0.1666 

0.2303 

0.2995 

0.3560 

0.5148 

0.5109 

-2 

0.001 1 

0.01 17 

0.0459 

o.i 101 

0.1656 

0.2287 

0.2971 

0.3554 

0 

0.0000 

0.0106 

0.0452 

0.1091 

0.1648 

0.2276 

0.2949 

0.3542 

■ 

0.5081 

2 

0.0000 ] 

0.0108 

' 0.0446 

0.1083; 

0.1640 

0.2269 

0.29.35 

] 0.3536 

] 0.5064 

4 

6 

o.oooo] 

o.oi io: 

[0.043 i] 

0.1067] 

0.1624 

0.2255] 

0.2919] 

03533 

0.5058 

o.oooo] 

0.0112 

0.041 1! 

0.1047; 

0.1604 

] 0.2240] 

0.2908 

6.3545 

] 0.5065 

8 

0.00 II 

0.0122 

0.0412; 

0.1045 

0.1603 

] 0.2243] 

0.2919] 

] 0.3594 

6.5085 

10 

0.0022; 

0 0124 

0.0419 

0.1043 

0.1605 

: 0.2238; 

0.2926j 

] 0.3652 

j 0.5122 

12 

0.0040 

i 

0.0127; 

0.0437 

i 

0.1042] 

0.1611 

] 0.2220] 

0.2924 

“63710: 

0.5183 

14 

0.0066 

0 0145 

0.0472] 

0.1046 

0.1623 

0 2217] 

0.2941 

0.3789] 

0.5267 

16 

0.0098 

0.0187 

0.0528] 

0.1058] 

0 1644 

1 0.2252' 

0.2998 

0.3908] 

6.5350 

18 

0.0136: 

0 0243 

0.0592; 

0.1083 

0.1678 

s 

0.2307 1 

1 

0.3080] 

0.4044; 

0 5408 

20 

0.0182; 

0.0302! 

0.0655 

0.1126, 

0.1727] 

0.2364] 

03166] 

0.4173] 

0.5447 

22 

0.0236; 

0.0362 

0.0716 

0.1182 

0.1788] 

0.2415: 

.1 

0.3254] 

0.4284] 

0.5476 

24 

0.0290 

0.0421 

0 0772 

0.1240 

6.1848 

0.2457’ 

i 

03337] 

0.4369 

0.5507 

26 

0.0337 

0.0472; 

0.0819 

0.1292 

0.1901] 

0.2499, 

0.3421] 

6.4443] 

0.5552 

28 

0.0379 

0.0514 

0.0857 

0.1340 

0.1948] 

0.2560 

0.3516 

0.4531 

0.5607 

30 

0.0421 

0.0549 

0.0888 

0.1383] 

0.1989] 

02652] 

0.3629; 

0 4581 

6.5608 

32 

0.0454 

0.0566 

0.0901 

0.1412] 

0.2017: 

0.2750- 

1 

63738; 

0.4665 

J 

0 5506 

34 

36 

0.047 1 

0.0565' 

0.0897 

0.1428* 

i 

0.2031 

1 

0.2836] 

0.3828, 

0.4821] 

0 5355 

0.0478 

0.0563 

0.0895 

0.1447] 

0.2049 

0.2910] 

0 3903; 

0.4829 j 

0 5221 

38 

40 

0.0482 

r- 

IT, 

o 

o 

0.0908 

0.1482] 

0.2085] 

02990] 

0.3986] 

0.4757] 

6.5109 

0.0489 

0.0598 

0.0934, 

0.1534] 

0.2138] 

0.3088! 

0.4084 

0.4583 

0.5007 
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Table 4. NACA pitot probe dynamic pressure calibration coefficient {C q ) values. 



131 , deg 

a , deg 

0 

5 

10 

15 

20 | 25 

30 

35 

40 

-40 

0.0440 

0.0473 

0.0489 

0.0604 

0.0366 

0.0147 

- 0.0159 

0.1368 

0.3529 

-38 

0.0129 

0.0118 

0.0148 

0.0249 

00118 

- 0.0108 

- 0.0764 

0.0299 

0.2597 

-36 

- 0.0183 

- 0.0201 

- 0.0160 

- 0.0068 

- 0.01 10 

- 0.0358 

- 0 . 1 03 ^ 

- 0.0473 

0.1799 

-34 

-32 

- 0.0477 

1 

O 

o 

\J\ 

o 

- 0.0475 

- 6.0391 

- 0.0381 

- 0.0587 

- 0.1242 

- 0.1038 

0.1032 

- 0.0705 

- 0.0755 

- 0.0749 

- 0.0676 

- 0.0663 

- 0.0728 

- 0.1432 

- 0.1308 

0.0300 

-30 

- 0.0835 

- 0.0911 

- 0.0946 

- 0.0894 

- 0.0907 

- 0.0857 

- 0.1302 

- 0.1485 

- 0.0351 

-28 

- 0.0937 

- 0.1034 

- 0.1 144 

- 0.1127 

- 0.1174 

. 

- 0.1143 

- 0.1280 

- 0.1614 

- 0.0936 

-26 

- 0.0911 

- 0.1021 

- 0.1242 

- 0.1269 

- 0.1356 

- 0.1387 

- 0.1372 

- 0.1806 

- 0.1424 

-24 

- 0.0851 

- 0.0966 

- 0.1219 

- 0.1299 

- 0.1430 

- 0.1525 

- 0.1462 

- 0.1894 

- 0.1776 



-22 

- 0.0804 

- 0.0919 

- 0.1133 

- 0.1289 

- 0.1466 

- 0.1615 

- 0.1560 

- 0.1948 

- 0.2043 

-20 

- 0.0690 

- 0 . 080.3 

- 0.1044 

- 0.1282 

- 0.1520 

- 0.1719 

- 0.1698 

- 0.1549 

- 0.2238 

-18 

- 0.0366 

- 0.0478 

- 0.0940 

- 0.1275 

- 0.1587 

- 0.1840 

- 0.1859 

- 0.1558 

- 0.2389 

-16 

- 0.0230 

- 0.0340 

- 0.0787 

- 0.1247 

- 0.1623 

- 0.1934 

- 0.1983 

- 0.1752 

- 0.2505 

-14 

- 0.0169 

- 0.0277 

- 0.0606 

- 0.1154 

- 0.1573 

- 0.1938 

- 0.2034 

- 0.1822 

- 0.2583 

-12 

- 0.0130 

- 0.0239 

- 0.0519 

- 0.1113 

-0 1557 

- 0.1968 

- 0.2149 

- 0.1949 

- 0.2674 

-10 

- 0.0101 

- 0.0207 

- 0.0468 

- 0.1082 

- 0.1546 

- 0.1989 

- 0.2262 

- 0.2106 

- 0.2773 

-8 

- 0.0071 

- 0.0177 

- 0.0447 

- 0.1070 

- 0.1551 

- 0.2015 

- 0.2356 

- 0.2263 

- 0.2866 

-6 

- 0.0048 

- 0.0152 

- 0.0446 

- 0.1070 

- 0.1568 

- 0.2044 

- 0.2428 

- 0.2401 

- 0.2947 

-4 

- 0.0026 

- 0.0129 

- 0.0445 

- 0.1068 

- 0.1578 

- 0.2063 

- 0.2471 

- 0.2502 

- 0.3014 

-2 

- 0.0008 

- 0.01 12 

- 0.0443 

- 0.1064 

- 0 . 1582 ; - 0.2074 

- 0.2495 

- 0.2566 

- 0.3069 

0 

0.0000 

- 0.0104 

- 0.0441 

- 0.1060 

- 0.1582 

- 0.2086 

- 0.2515 

- 0.2612 

- 0.3117 

2 

0.0000 

- 0.0106 

0.0436 

i J 

- 0. 1053 

- 0.1576 

- 0.2094 

- 0 . 25.32 

: - 0.2647 

- 0.3154 

4 

0.0000 

- 0.0108 

- 0.0421 

- 0.1036 

- 0.1 559 

- 0.2092 

- 0 . 25391 - 0.2671 

: -03 184 

6 

0.0000 

- 0.0110 

- 0.0402 

- 0.1016 

- 0.1538 

- 0.2086 

- 0.2541 

- 0.2695 

1 - 0.3208 

8 

- 0.0011 

- 0.0120 

4 ) 3)402 

-0 1014 

- 0.1537 

- 0.2093 



- 0.2556 

- 0.2741 

' - 0.3227 

10 

12 

- 0.0021 

- 0.0122 

- 0.0409 

- 0.1012 

- 0.1538 

- 0.2087 

- 0.2558 

- 0.2780 

- 0.3247 

- 6.0038 

- 0.0124 

- 0.0426 

- 0. 1011 

- 0.1542 

- 0.2062 

- 0.2539 

- 0.2802 

- 0.3271 

14 

- 0.0063 

- 0.0140 

- 0 . 0460 : 

- 0. 1014 

- 0.1551 

- 0.2046 

- 0.2528 

- 0.2828 

; - 0.3299 

16 

- 0.0092 

- 0.0180 

- 0 . 0512 : 

- 0.1022 

- 0.1564 

- 0.2061 

- 0.2548 

- 0.2877 

U 0.3307 

18 

- 0.0127 

s - 0 . 0233 : 

- 0.0573 

- 0.1040 

- 0.1586 

. 

- 0.2093 

- 0.2582 

.- 0.2928 

| - 0.3273 

20 

- 0.0170 

i - 0.0288 

- 0.0630 

0.1072 

- 0.1617 

- 0.2119 

- 0.2607 

- 0.2955 

1 - 0.3201 

22 

- 0.0220 

; - 0.0344 

- 0.0684 

i - 0 . 1 114 

- 0.1653 

- 0.2129 

- 0.2617 

- 0.2947 

1 - 0.3100 

24 

26 

28 

- 0.0268 

| - 0.0396 

- 0.0730 

l - o.l 152 

- 0.1678 

- 0.2116 

. 

- 0.2601 



- 0.2896 

| - 0.2982 

- 0.0305 

. - 0.0437 

- 0.0762 

; - 0.1 178 

- 0.1686 

- 0.2083 

- 0.2561 

- 0.2814 

f - 0.2856 

- 0.0332 

- 0.0464 

- 0.0778 

-0 1191 

- 0.1674 

- 0.2037 

- 0.2499 

- 0.2717 

- 0.2710 

30 

- 0.0349 

- 0.0473 

- 0.0773 

1 - 0.1186 

- 0.1636 

- 0.1981 

- 0.2415 

- 0.2546 

- 0.2472 

32 

- 0.0339 

- 0.0444 

- 0.0729 

- 6.1 142 

- 0.1551 

- 0.1899 

- 0.2299 

- 0.2381 

- 0.2101 

34 

- 0.0298 

\ 

o 

o 

00 

o 

- 0.0650 

- 0.1048 

- 0.1404 

- 0.1778 

L.„ - — „ ...„ 

- 0.2144 

- 0.2269 

- 0.1661 

36 

- 0.0224 

- 0.0289 

- 0.0544 

- 0.0925 

- 0.1217 

- 0.1619 

- 0.1959 

- 0.1992 

- 0.1221 

38 

- 0.0104 

- 0.0161 

- 0.0402 

- 0.0772 

- 0.1016 

- 0.1441 

- 0.1768 

- 0.1622 

- 0.0788 

40 

0.0086 

0.0030 

- 0.0198 

- 0.0561 

- 0.0775 

- 0.1222 

- 0.1541 

- 0.1098 

- 0.0315 
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Table 5. Airplane axis coordinates of system components. 


A, ft 

3*. 0 

2, ft 

Description 

-6.3675 

19.7242 i 

-2.6875 

right boom 

-6.3658 ‘ 
5.6108 ' 

-19.7242 i 
2.1308 

-2.6875 

0.6292 

left boom tip 
nose boom tip 

14.5833 [ 

1.2083 i 

0.5779 

body axis origin (a^, y^, z^) 


Table 6. OVDRA calculated and recorded parameters. 


Symbol 

Units 


Eqn # 

Description 

a r 

deg 


22 

; Right boom local angle of attack 

«/ 

deg 

\ 

23 

: Left boom local angle of attack 

a n 

deg 

1 

24 

: Nose boom local angle of attack 

Or 

deg 


25 

Right boom flank angle 

0/ 

deg 

\ 

26 

Left boom flank angle 

Pn 

deg 

j 

28 

Nose boom local sideslip angle 

Qr 

psi 

| 

38 

Right boom dynamic pressure 

Qi 

psi 

L. 

38 

Left boom dynamic pressure 

Qn 

psi 


39 

Nose boom dynamic pressure 

Ps r 

psi 


36 

Right boom static pressure 

Psj 

psi 

— h 

36 

! Left boom static pressure 

Ps„ 

psi 


37 

Nose boom static pressure 

M r 

- 

\ 

40 

Right boom local Mach number 

M f 

- 

i 

40 

i Left boom local Mach number 

M n 

- 


40 

Nose boom local Mach number 

T abs . 

°R 


42 

Local absolute temperature 


n R 


43 

Local free-stream air temperature 

T e 


\ 

44 

Potential temperature 

a 

ft/sec 

| 

45 

Local speed of sound 

V*loc r 

ft/sec 

■ ’i"’ 

46 

Right boom local true airspeed 

Vtyocf 

ft/sec 

\ 

46 

Left boom local true airspeed 

Vt loc„ j 

ft/sec 

r 

? 

i 

46 

Nose boom local true airspeed 

Pr 

deg 

■* 

27 

Right boom local sideslip angle 

Jl 

deg 

1 

21 

Left boom local sideslip angle 

l % 1 

ft/sec 

""T 

I 

! 

i 

47 

Right boom body axis longitudinal airspeed 

«<v j 

ft/sec 

1 

j... 

47 

Left boom body axis longitudinal airspeed 

Uh 

u n } 

ft/sec 

1 

47 

: Nose boom body axis longitudinal airspeed 
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Table 6. Continued. 


Symbol 

Units 

Eqn # 

Description 


ft/sec 

48 

Right boom body axis lateral airspeed 

V/,, 

ft/sec 



48 

Left boom body axis lateral airspeed 

x \ 

ft/sec 

48 

Nose boom body axis lateral airspeed 

n v 

ft/sec 

49 

Right boom body axis vertical airspeed 


ft/sec 

49 

Left boom body axis vertical airspeed 

M \ 

ft/sec 

49 

Nose boom body axis vertical airspeed 

Vt r 

ft/sec 

50 

Right boom true airspeed 

Vt, 

ft/sec 

50 



Left boom true airspeed 

Vt n 

ft/sec 

50 

Nose boom true airspeed 

a 'r 

deg 



51 

Right boom angle of attack 


deg 



51 

Left boom angle of attack 

a t 

n 

deg 

51 

Nose boom angle of attack 

P'r 

deg 

52 

Right boom sideslip angle 

P‘, 

deg 

52 

Left boom sideslip angle 

P‘,, 

deg 

52 

Nose boom sideslip angle 

Vt avg 

ft/s ec 

53 

Average true airspeed 

“V s \ 

ft/sec 

54 

Average body axis longitudinal airspeed 

l 'W i 

ft/sec j 

i 

: 55 

Average body axis lateral airspeed 

H’Z, 

ft/sec j 

56 

Average body axis vertical airspeed 


deg 

57 

Average angle of attack 

Pavg 

deg 

58 

Average sideslip angle 

T*(U) 

1 

| 

60 ; 

Body to Earth axis transformation matrix element 

T^d.2) ; 

| 

60 

Body to Earth axis transformation matrix element 

T*(l,3) ; 

| 

60 

Body to Earth axis transformation matrix element 

"•'/)<> (2. 1 ) 

! 

60 

Body to Earth axis transformation matrix element 

^ (2,2) 

” j 

60 ; 

Body to Earth axis transformation matrix element 

(W) j 

1 

| 

60 1 

Body to Earth axis transformation matrix element 

Tie (3,1) ! 

! 

60 

Body to Earth axis transformation matrix element 

Tie (3.2) I 

! 

60 

Body to Earth axis transformation matrix element 

T*<3.3) ; 

\ 

60 | 

Body to Earth axis transformation matrix element 

W Nr \ 

ft/sec | 

I 

61 

Right boom northerly wind velocity 

W Er 1 

ft/sec j 

62 ! 

Right boom easterly wind velocity 

W'Zr : 

ft/sec j 

63 

Right boom vertical wind velocity 
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Table 6. Concluded. 


Symbol 

Units 


Description 

W N , 

ft/sec 

61 

Left boom northerly wind velocity 

w E i 

ft/sec 

62 ! Left boom easterly wind velocity 

w Z i 

ft/sec 

63 

Left boom vertical wind velocity 

w Nn 

ft/sec 

61 

' - 

Nose boom northerly wind velocity 

W E n 

ft/sec 

62 Nose boom easterly wind velocity 

W Z n 

ft/sec 

63 Nose boom vertical wind velocity 

Wj\[ avg 

ft/sec 

61 Average northerly wind velocity 

W E civ g 

ft/sec 

62 Average easterly wind velocity 



W'Zavg 

fl/sec 

63 Average vertical wind velocity 

■ ..... .... . . r _ . . — 

V wnd 

knots 

66 Average wind speed 

Vwnd i de S 

I 67 Average wind direction in degrees true 






Wake & Atmospheric 
Measurement 



Figure 1 . - Wake measurement flight test setup. 



Figure 2. - C-130 wake generator airplane. 
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Figure 3. - OV-IO major instrumentation systems 



Figure 4. - OV-IO tail video of wake measurement. 
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(a) Pallet outside of cargo bay. (b) Pallet installed for flight. 

Figure 5. - Data acquisition system research pallet. 



Figure 6. - Research equipment in the forward cargo-bay area. 
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Figure 7. - Right boom NACA pitot probe with angle of attack and sideslip vanes. 



Figure 8. - Nose boom 5-hole probe. 
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Latitude 


Figure 11.- Euier angles with body and inertial axis systems. 




Check which vortex encountered?^ 


Vrtx = R no 


Left vortex 
^=^CI30-90 


Right vortex 
'V/ip ='* , C130 +9 ° 



Compute location of wing tip: 

Lett tip , Lng up = f{Lat C \ 30 (')> Lng Cl 3o(<). ^tip < \ h C\ 30 ) 


Compute location of drifted wake: 

Lat w k , Lng w f, — f \l-Mt , Lng tl ^ , 'P Vi7W / ~ 1 80 , d w f, ) 
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n 



46 









Appendix A - OV-IO Airspeed Calibration 


The objective of the OV-IO airspeed calibration was to determine the position-error pressure 
corrections and the angle-of-attack and sideslip calibration for the three booms. The airspeed calibration 
flights were Flight 613 flown on September 24, 1996 and Flight 614 flown on November 13, 1996. The 
position error corrections for the left, right and nose boom were computed in a similar manner as the 
“camera method” described in reference 5. The difference between the method outlined here and the 
camera method was that the airplane altitude was determined from differential GPS measurements rather 
than photographically. 

Flight Test Procedure 

The basic test procedure involved flying a series of steady, low altitude, level passes at different 
airspeeds and sideslip angles. A row of 140-foot tall light towers at Langley Airforce Base (LAFB) 
served as a visual reference for the level passes. After every two or three low altitude passes the OV-IO 
would land and collect static surface data such as temperature, pressure and wind measurements. Pre- and 
post-flight instrument calibration data and in-flight steady trim data were also collected to determine 
instrumentation bias values. The data from the OV-IO experimental data system was reduced from 
128 Hz to I Hz using a one-second-averaging interval. 

Computing Reference Values 

The first step in determining the calibration corrections is to compute the reference values from which 
the measured values are calibrated. This basically involves extrapolating surface measurements to the 
pass altitude based on a measured lapse rate. 

Reference Angle of Attack and Sideslip Calculation 

The wind velocity components during the pass must be determined to compute the reference sideslip 
and angle of attack. The LAFB surface winds, which were measured 12 feet above the runway, were 
recorded throughout each flight. Figure A I shows the surface wind time history for both flights. The 
numerous takeoffs and landings throughout the flight provided a time history of the wind change with 
altitude. The change in the wind speed and direction with altitude was computed for each takeoff and 
landing through a linear fit of the measured winds from 8 to 16 feet and from 120 to 160 feet above the 
runway. The assumption here is that although the magnitude of OV-IO wind velocity and direction 
measurements were not yet calibrated and therefore not valid, the measured change of the wind velocity 
and direction were essentially correct. The time histories of the computed wind gradients for both flights 
are shown in figure A2. 

The reference wind speed and direction for a given pass was computed through a linear interpolation 
of the surface wind speed and direction and the altitude-gradients as a function of the pass time and 
altitude. The altitudes were determined from the differential Ashtech GPS measurements. 

^wnd pass ' = fp pass ' ' pass ~ ^ grnd ) + ^wndg rn j ^ pass ) (knots) (A 1 ) 

^wnd^pass^^fT-^pass^pass h grnd) + ' ¥ W nd grnd « pass'* (deg- magnetic) (A2) 
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The inertial wind components are computed from the wind direction and speed through: 


W N =-\M8 V wnJ cos| 


W £ =-1.688V^sin 




wnd. 


mag 


* U.'l 


wruJ„ 


mag 

W z =0 (fps) 

where the magnetic variation ( Mag vtir ) at the test site was 9.23 degrees west. 
The reference true airspeed was computed as: 


- Ma Svar 

) 

| (fps) 

(A3) 

\ 

— Mag var 

/ 

(fps) 

(A4) 



(A5) 


Vt ref = iI(n -wJ+(e- w e f +(-Z-W z f (fps) 


The air-mass-relative body-axis velocity components were computed from: 


(fps) 


f 



'V 


n-w n ' 

V 

it 

e-w e 

M V 

V J ) 


1 

N 

1 


where T* is defined in equation 60. 

The reference angle of attack was computed as: 


°V = tan 


-i 


f w„ ' 


U- 


b,. 4 


(deg) 


and the reference sideslip from: 


P, e/ =sin 


-1 


ref J 


/V ^ 

b ref 


Vt 


ref 


(deg) 


(A6) 


(A7) 


(A8) 


(A9) 


The wing boom balsa vanes measure the flank angle rather than the sideslip angle. The reference 
flank angle was computed as: 


tan 


"ref 


■ tan 


P ref 


cos a 


(deg) 


(A 10) 


ref 
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Reference Static Pressure Calculation 

The reference static pressure was determined by extrapolating the surface measurements to the probe 
altitude assuming a standard density lapse rate. Time histories of the ground static pressure and 
temperature measurements from the OV-IO are shown in figures A3 and A4, respectively. Also shown in 
figure A3 are the static pressure measurements from an independent ground measurement. The altitude of 
each probe was determined from the differential Ashtech GPS measurement and the INS attitude 
measurement. 

hj = 3.281 * HI -\x b -xj jsin 0 93 - \j jcos 0 93 sin $93 + -Zj )cos0 93 cos<t) 93 (ft) 

(All) 

j -l, n 9 r 


The reference static pressure was computed from: 

Ps ref,j = Ps grnd ~ 8 ' Pnud.j ( /; j ~ h grnd )/ 144 <P si ) 
j = in, r 


(A 12) 


The density at the mid-altitude point was computed as: 

P midj ~ P grnd + "2 ~ ^ grnd ) 

j = l, n 9 r 


(A 13) 


with — = -6.8x10 7 slugs/ft 4 (standard density lapse rate). The air density at the ground was corrected 
dh 


for non-standard conditions through: 
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with /? = 1716.5 
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Compute Position Error Correction 

The position error correction was determined via a least squares linear-fit of the difference between 
the probe and reference static pressures as a function of the probe dynamic pressure. 

Probe Static and Dynamic Pressure Calculation 

The probe static and dynamic pressures were corrected for local flow angularity based on wind tunnel 
derived calibration coefficients. The static and dynamic pressure coefficients are defined as: 

c s = PSref PSprohe (A 15) 

Q ref 
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and 


c = V q probe (A|6) 

“ 4 ref 

The reference values (subscripted ref) were the wind tunnel “truth” measurements. For the wing tip 
mounted NACA probes the coefficients were determined via a two-dimensional table lookup as a 
function of the reference angle of attack and sideslip. Tables 3 and 4 list the C s and C (j table lookup 

values. The corrected wing tip probe dynamic (q c ) and static ( Ps t . ) pressures were computed from 
equations 29 and 30, respectively. The nose boom 5-hole probe dynamic ( q 5 ) and static (Ps 5 ) pressures 
were computed from equations 20 and 21, respectively. 

The position-error pressure correction was computed as the difference between the calibrated probe 
measurement and the reference value. 


A Pq j = Ps c j - Ps re f j (psi) 

j = lr 

AP ‘l n = 5 — ref ,n 

Figures A5, A6, and A7 show the position error values of the nose, left and right boom probes, 
respectively, as a function of the probe calibrated dynamic pressure for the zero sideslip passes of both 
flights. Also shown on the figures is the least squares linear fit to the data. The linear fit coefficients 
were used as the position error correction coefficients for zero sideslip. 

The effect of sideslip on the position error was determined by subtracting the zero sideslip position 
correction from the position error values at nonzero sideslip. 

A Pq p = A Pc/ ( p^oj - A Pq f „( p=o) (A 1 9) 

Figures A8, A9, and A 10 show the position error values due to sideslip of the nose, left and right boom 
probes, respectively. Also shown on the figures is the least squares linear fit to the data. The linear fit 
coefficients were used as the position error correction coefficients for non-zero sideslip. For the nose 
boom probe the effect of sideslip was symmetric about zero sideslip. For the wing tip probes the linear fit 
varied for positive and negative sideslip angles. The left and right boom corrections were nearly mirror 
images of each other. Figure A 1 1 shows the two data sets combined with the left being the opposite sign 
as the right. The linear fit of the combined data sets was used for the tip boom correction coefficients. 

The position error correction is applied to the static and dynamic pressure measurements through 
equations 36 through 39. 

Angle of Attack Calibration 

The angle of attack calibration was determined via a least squares linear fit of the measured value as a 
function of the reference value obtained from equation A8. The nose boom measured angle of attack was 
determined from equation 18. The wing tip values were obtained from the balsa vane measurements. 


(A 17) 
(A 18) 
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The airplane angular rates were assumed to be negligible. The data and linear fits are shown in figures 
A 12, A 13, and A 14 for the nose, left and right booms, respectively. There was an apparent offset in the 
wing tip angle of attack measurements between flights 613 and 614. This was presumably due to a slight 
warping of the balsa vane with variation in temperature and humidity. Due to this offset the slope of the 
right and left boom calibration was computed from flight 613 values only. Flight 613 had a larger 
variation in the range of angles than Flight 614. The wing tip offset values for flight 614 were then 
computed using the slope calculations from flight 613. 

Sideslip Angle Calibration 

The sideslip calibration was determined in a similar manner as the angle of attack calibration. The 
reference sideslip value was obtained from equation A9. The nose boom measured sideslip angle was 
determined from equation 19. Figure A15 shows the nose boom sideslip calibration. There was a slight 
unresolved offset in the nose boom measurements between flights 613 and 614. The slope of the nose 
boom calibration curve was determined from flight 614 data since the majority of the sideslip data was 
obtained on that flight. 

Figures A16 and A 1 7 show the flank angle data and calibration for the left and right boom, respectively. 
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Figure A4. Ground temperature history for flights 613 and 614. 
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Figure A5. Nose boom 5-hole probe position error correction. 
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Figure A6. Left wing tip boom probe position error correction. 
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Figure A7. Right wing tip boom probe position error correction 
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Figure A8. Nose boom probe position error correction with sideslip. 
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Figure A9. Left boom probe position error correction with sideslip. 
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Figure A10. Right boom probe position error correction with sideslip. 
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Figure A1 1 . Combined right and left probe position error correction with sideslip 
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